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Abstract 

We discuss and compare different approaches to include gluon 
transverse momenta for heavy quark-antiquark pair production. The 
correlations in azimuthal angle and in heavy quark, heavy antiquark 
transverse momenta are studied in detail. The results are illustrated 
with the help of different unintegrated gluon distribution functons 
(UGDF) from the literature. We compare results obtained with on- 
shell and off-shell matrix elements and kinematics and quantify where 
these effects are negligible and where they are essential. We concen¬ 
trate on the region of asymmetric transverse momenta of charm quark 
and charm antiquark. Most of UGDFs lead in this corner of the phase 
space to almost full decorellation in azimuthal angle. We propose 
correlation observables to be best suited in order to test the existing 
models of UGDFs. 

PACS: 12.38.-t, 12.38.Cy, 14.65.Dw 


1 Introduction 

The charm-anticharm production in photo- and hadroproduction is known 
as one of the crucial tests of conventional gluon distributions within a stan¬ 
dard factorization approach. At high energies one tests gluon distributions at 
low values of longitudinal momentum fraction. At leading order of collinear 
approach the heavy quark and heavy antiquark are produced back to back 
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i.e. the azimuthal correlation function is proportional to S{(j) — tt) and the 
distribution in pi ^ (heavy quark transverse momentum) and p 2 ^t (heavy an¬ 
tiquark transverse momentum) is proportional to — P 2 ,t)- In collinear 
approach this changes only in next-to-leading order. In unintegrated gluon 
distribution approach the azimuthal angle and pt decorrelations are obtained 
already in leading order of perturbative expansion. 

Different models of UGDF were presented in the literature. Recently we 
have tested some of the models in photoproduction [T]. Here we shall discuss 
in detail the case of hadroproduction. 

It was realized recently that the understanding of transverse momenta can 
be a clue to understand many high energy reactions. In some reactions like 
meson hadroproduction, prompt photon and Drell-Yan production or single 
spin asymmetry (SSA) the transverse momenta are not obtained dynamically 
but are included “by hand” with Gaussian smearing (see e.g. IZIE])- The 
parameter of the smearing is assumed to be independent of longitudinal mo¬ 
mentum fraction and is usually arbitrarily chosen or adjusted to a selected 
set of data. 

In some of high-energy reactions that involve gluons the unintegrated 
gluon distributions are modelled or are obtained as a solution of some (usually 
simplihed) QGD evolution equations. In particular, the charm-anticharm 
hadroproduction is considered as one of the flag processes for the so-called 
/c^-factorization approach lUEiin]. There are already a few phenomenological 
trials to calculate the charm or bottom production within the /c^-factorization 
approach 0 012!. The results of these trials differ considerably due to 
completely different UGDFs used. 

In the present paper we shall compare different approaches how to in¬ 
clude transverse momenta proposed in the literature for the charm-anticharm 
pair production at the Tevatron energy W = 1.96 GeV both for the single 
charm/anticharm distributions as well as for a few charm-anticharm correla¬ 
tion observables. 


2 Formalism 

Let us consider the reaction hi + h 2 —>■ Q + Q + X, where Q and Q are 
heavy quark and heavy antiquark, respectively. In the leading-order (LO) 
approximation within collinear approach the triple-differential cross section 
in rapidity of Q (yi), in rapidity of Q ( 2 / 2 ) and transverse momentum of one 
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of them [pt) can be written as 


da 

dyidy2d‘^Pt 


1 


'^xipi{xi,y^) X2Pj{x2,y^) 

ij 


( 1 ) 


Above Pi(a:i,/i^) and Pj{x 2 ,p‘^) are familiar (integrated) parton distribntions 
in hadron hi and ^ 2 , respectively. There are two types of LO 2 —> 2 
snbprocesses which enter Eq. (P): 99 QQ cind qq —> QQ. The hrst 
mechanism dominates at large energies and the second near the thresh¬ 
old. The parton distribntions are evalnated at: Xi = ^ (exp(?/i) -|- exp(?/ 2 )), 

X 2 = ^ (exp(—yi) -|-exp(— 1 / 2 )), where rrit = \Jp^ + The formnlae for 
matrix element sqnared averaged over initial and snmmed over hnal spin 
polarizations can be fonnd e.g. in Ref. PI- 

If one allows for transverse momenta of the initial partons, the snm of 
transverse momenta of the final Q and Q no longer cancels. Formnla o can 
be easily generalized if one allows for the initial parton transverse momenta. 
Then 


da 

dyidy2d‘^Pi^t(Pp2,t 



d'^Ki^t (PK2,t 1 
71 71 167l^{XiX2Sy 


\A4.ij 
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(5^ + /? 2 ,t - Pl,t - P 2 ,t) Ti{xi, k{ ^) Tj{x 2 , ni t) , 


( 2 ) 


where now Ti{xi, k{ J and J^j{x 2 , ^ 2 1 ) so-called nnintegrated parton dis¬ 
tribntions The extra integration is over transverse momenta of the initial 
partons. The two extra factors I/tt attached to the integration over d^Ki^t and 
d‘^K 2 ^t instead over dn^ ^ and dK^ ^ as in the conventional relation between nn¬ 
integrated (JF) and integrated {g) parton distribntions. The two-dimensional 
Dirac delta fnnction assnres momentnm conservation. Now the nnintegrated 
parton distribntions mnst be evalnated at: xi = (exp(?/i)-|-exp(|/ 2 )), 

X 2 = ^ (exp(-|/i) exp(-|/ 2 )), where rrii^t = \JpXt + general, the 

matrix element mnst be calcnlated for initial off-shell partons. The corre¬ 
sponding formnlae for initial glnons were calcnlated in Pd (see also [HI)- In 
the present paper for illnstration we shall compare resnlts obtained for both 
on-shell and off-shell matrix elements. It is easy to check that in the limit 

^In this paper we shall use parallel two different conventions of unintegrated gluon 
distribjations: 

(a) fg JF(x, K^)dK^ xg(x, y,^), 

(b) fo^ f(x,K^)^ ~ xg{x,y^). 
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where t, u, s are Mendelstam variables for the partonic subprocess. 
Introducing new variables: 

Qt = + i^2,t , 

qt = «i,t - « 2 ,t 


( 4 ) 


we can write: 


dcTij 


dyidy2(Ppi,td?P2,i 


= / d'^qt 


IM 


47r2 167r2(a;ia;2s)2 
Ti{xi,Kly) J^j{x 2 ,K\t) ■ 


( 5 ) 


This formula is very useful to study correlations between the produced heavy 
quark Q and heavy antiquark Q. 

For example 

d(Jij 

dpi,tdp2,t 


In the last equation we have introduced 0. = 01 — 02, where 0_ G (-27r, 
27r). The factor 4 tt comes from the integration over 0_,_ = 0i + 02. The hrst 
factor 1/2 comes from the Jacobian transformation while the second factor 
1/2 takes into acconnt an extra extension of the domain when using 0+ and 
0 _ instead of 0i and 02 . 


d0id02 pi,tP2,t / dyidy2 / d^qt 


J^i{xi,Kly J^j{x2,K,lt) 

47r ^ i y d(l)_ pi,tP2,t j dyidy2 J d‘^qt 

J^i{xi,Kly J^j{x2,Kl t) . 


47r2 167r2(a;ia;2s)2 

1 1 


\M, 


\M, 


47r2 167r2(a;ia;2s)2 

( 6 ) 
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At the Tevatron energy the contribution of the gg —>■ QQ subrocess is 
more than order of magnitude larger than its counterpart for the qq —*■ QQ 
subprocess. Therefore in the following we shall take into account only gluon- 
gluon fusion process i.e. i=0 and j=0. The quark-antiquark mechanism is 
important only at small center-of-mass energies. 

In LO collinear approach 


da 

dpi,tdp2,t 


ot d{pi,t-P2,t) ■ 


(7) 


A deviation from this relation is therefore a measure of the initial gluon 
(parton) transverse momenta. 

Purely perturbative ^ /crfactorization formalism to hih 2 —>■ QQ applies 
if Kit,>^ 2 t > ^ 0 - The choice of Kq is to large extent arbitrary. In Refs.[7j 
a rather large Kq was chosen and the space x ^ was subdivided into 
four disjoint regions. For example the contribution when both nf ^ and t 
are small was replaced by the leading-order collinear cross section. Such an 
approach by construction assures that a^Q > (collinear LO). 

It is rather obvious that the total cross section strongly depends on the 
choice of Kq. Our philosophy here is different. Many models of UGDF in the 
literature treat the soft region explicitly. Therefore we use the /c^-factorization 
formula everywhere on the x plane. For perturbatively generated 
UGDFs, like KMR ^21 for instance, we shall extrapolate the perturbative 
behaviour into the soft region. Thus our approach guarantees a smooth 
behaviour on the x plane and reduces arbitrariness in dividing the 
two-dimensional gluon-transverse-momentum space. 


3 Unintegrated gluon distributions 

Obtaining UGDFs from underlying QGD is not an easy task. The main rea¬ 
son is the difficulty in separating out the perturbative from nonperturbative 
domains. The nonperturbative domain is difficult per se. Even in pertur¬ 
bative region different schemes of resummation have been proposed. In this 
section we describe briefly some representative UGDFs from the literature 
used in the present paper to calculate charm-anticharm production. 

3.1 Kwiecinski gluon distribution 

Kwiecihski has shown that the evolution equations for unintegrated parton 
distributions takes a particularly simple form in the variable conjugated to 

^when both UGDFs are generated perturbatively 


5 




the parton transverse momentum. The two possible representations are in¬ 
terrelated via Fourier-Bessel transform 


Tk{x,K;,ij,p) = / dbbJo{Ktb)J^k{x,b,iJ.p), 


Tk{x,b,llp)= / dKt KtJo{Ktb)J^k{x,K^,IJ,p) . 


( 8 ) 


The index k above numerates either gluons (k=0), quarks (k> 0) or anti¬ 
quarks (k< 0). 

In the impact-parameter space the Kwiecihski equation takes a rather 
simple (diagonal in b) form |TT1[T21II2] (see also [TK]b 

The perturbative solutions b, /i|.) do not include nonperturbative 

effects such as, for instance, intrinsic transverse momenta of partons in collid¬ 
ing hadrons. One of the reasons is e.g. internal motion of constituents of the 
proton. In order to include such effects we modify the perturbative solution 
j^pert(^, write the modihed parton distributions ^g{x, b, /i|,) in the 

simple factorized form 


b, A) = b, mP ■ F^^ib) . (9) 


In the present study we shall use the following functional form for the non¬ 
perturbative form factor 

F7((>) = = exp . (10) 

In Eq. (HHD&o is the only free parameter. 

In the following we use leading-order integrated parton distributions from 
Ref.j2Sj as the initial condition for QCD evolution. The set of integro- 
differential equations in b-space is solved by the method based on the dis¬ 
cretisation made with the help of the Chebyshev polynomials (see e.g. na). 
Then the unintegrated parton distributions jFg(x, 6,are put on a grid in 
X, b and /i^ and the grid was used in practical applications for Chebyshev 
interpolation. In our practical application we need rather gluon distributions 
in momentum space. The latter are obtained via Fourier transform (jH)) from 
Tg{x,b,fJ,l^). 


3.2 Kimber-Martin-Ryskin distribution 

Resumming virtual contributions to DGLAP equation, the unintegrated par- 
ton distributions can be written as UTI 

fa(x, lip) ^ Ta{K^, pp ■ ^ ^ f (j) (^) o' (j <iz . (11) 

n' ^ 
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Specializing to the gluon distribution the Sudakov form factor reads as 


= exp - 


dp^ as {p^ 


- 1-5 


^2 p^ 271 


dz 




( 12 ) 

The Sudakov form factor introduces a dependence on a second scale 
Different prescriptions for 6 have been used in the literature. The most 
popular choice 6 = pf/{pf + i^) corresponds to the angular ordering in the 
gluon emission. 

It is reasonable to assume that the unintegrated gluon density given by 
(HH) starts only at |THI, i.e. in the perturbative domain. At lower 

nf an extrapolation is needed. A use of the GRV integrated gluon distribu¬ 
tion 1211 EE] in m seems more adequate than any other gluon distribution 
because it allows to go down to small gluon transverse momenta. Following 
Ref. M = 0-5 GeV^ is taken as the lowest value where the unintegrated 
gluon distribution is calculated from Eq. (ED. Below it is assumed 


J^{x, = f{x, = f{x, nl)/nl . (13) 


In general, the quantity k^q can be treated as a free nonperturbative param¬ 
eter. 

The choice of /i|, in our case of charmed quark-antiquark production is not 
completely obvious. In the present analysis, because we limit to not too large 
quark/antiquark transverse momenta, pjp = is taken for simplicity. 

This allows to prepare a two-dimensional grid J^KMR{xi, k|) to be used for 
further interpolation in the /ci-factorization formula which accelerates the 
multi-dimensional integrations. 

If Tg in Ea. (lll|l is ignored we shall denote the corresponding gluon dis¬ 
tribution as foGLAP or Tdglap and call it DGLAP gluon distribution for 
brevity. 


3.3 BFKL gluon distribution 

At very low x the unintegrated gluon distributions are believed to fulhl BFKL 
equation HSj. After some simplihcations I2DI the BFKL equation reads 

df{x,q^) ^ asNc 2 f{x,qu) " f{x,qt) 

^ dx 71 Jo qlt [ \q^-qu\ ^/qt + Mlt 

(14) 

The homogeneous BFKL equation can be solved numerically m- Here in the 
practical applications we shall use a simple parametrization for the solution 
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EH 

fix, «?) 


x^\qiJ ^27rA"ln(l/a:) [ r2A"ln(l/a;)_ 


(15) 


In the above expression A = 4Q;sln2, \" = 28 a<j(C(3), = Stts/vr, ^(3) = 

1.202. The remaining parameters were adjusted in HQ to reproduce with 
a satisfactory accuracy the gluon distribution which was obtained in HQ as 
the numerical solution of the BFKL equation. It was found that q = Qq = 1, 
C(j)o = 1.19 and r = 0.15 EH- 


3.4 Golec-Biernat—Wiisthoff gluon distribution 

Another parametrization of gluon distribution in the proton can be obtained 
based on the Golec-Biernat-Wiisthoff parametrization of the dipole-nucleon 
cross section with parameters htted to the HERA data. The dipole-nucleon 
cross section can be transformed to corresponding unintegrated gluon distri¬ 
bution. The resulting gluon distribution reads EH- 

asJ^ix, expi-Rlix)Kt) , (16) 

where 

- ifT (fv ■ 

From their £t to the data: cxo = 29.12 mb, xq = 0.41 ■ 10“^, A = 0.277. In 
order to determine the gluon distribution we take somewhat arbitrarily = 

0 . 2 . 


3.5 Gluon distribution a la Kharzeev-Levin 


Another parametrization, also based on the idea of gluon saturation, was 
proposed in |23|. In contrast to the GBW approach ra. where the dipole- 
nucleon cross section is parametrized, in the Karzeev-Levin approach it is 
the unintegrated gluon distribution which is parametrized. In the following 
we shall consider the most simplihed functional form: 


Rix, K^) 


/o 

/o ■ I? if > Q'l 


(18) 


The saturation momentum Qg is parametrized exactly as in the GBW model 

= 1 GeV2 . (a)\ 

The normalization constant /o was adjusted in EE! to roughly describe 
the HERA data. 
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3.6 Naive Gaussian smearing 

For better understanding we shall compare our results with the results ob¬ 
tained with a simple unintegrated gluon distribution: 

Tnaive{^-, /i|) ■ fcaussi^^) , (19) 

where is a standard collinear (integrated) gluon distribution and 

fcaussii^'^) is a Gaussian two-dimensional function: 

fcaussiK^) = {-'^V‘^^1) A • (20) 

Such a phenomenological procedure is often used to improve collinear calcu¬ 
lations for small transverse momenta j2]. The UGDF dehned by Ea. dlbj) and 
(PU|) is normalized such that: 

j ^naive{x, /i^) dti^ = Xg''°^\x, /i|) . (21) 

4 Results 

In the present analysis we shall discuss both inclusive spectra of charm 
quarks/antiquarks as well as correlations between quark and antiquark. In 
order to obtain the integrated cross section a 7-dimensional integration has 
to be performed. We have carefully explored the 7-dimensional phase space 
in order to optimize the integration. If not otherwise stated we take the 
whole range of quark/antiquark rapidities: -6 < |/i < 6 and -6 < ?/2 < 6. 

Let us start from single particle spectra of charm quarks. In Fig^ we 
demonstrate the effect of gluon transverse momentum Gaussian smearing on 
hnal quark transverse momentnm distribution for W = 1.96 TeV, i.e. at 
the present Tevatron energy. The smearing of primordial gluon distributions 
causses only a mild broadening of the charm quark-antiquark spectra. There 
is almost no difference if off-shell kinematics is used instead of on-shell one 
that is usually used in the case of the ad hoc Gaussian smearing. 

In FigOwe show a similar effect for the Kwiecihski UGDF for hxed fac¬ 
torization scale /i|^ = 4m^. There is rather weak dependence on the value of 
6o used for the nonperturbative form factor. The dependence on factorization 
scale is shown in Figl^l The resulting effect is rather mild. 

The off-shell matrix elements necessary for the /c^-factorization approach 
were calculated only for some selected reactions. According to our knowledge 
it was never shown quantitatively in the literature what is actual difference if 
the off-shell matrix element is replaced by its on-shell counterpart. In Fig 13] 
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we present results for the Kwiecihski UGDF for both on-shell and off-shell 
matrix elements. For this observable the use of the off-shell matrix element 
causes only a small enhancement compared to the use of the on-shell matrix 
element. 

In FigElwe collected results obtained with different unintegrated gluon 
distributions from the literature. In this case consequently the off-shell ma¬ 
trix element and off-shell kinematics were used. The GBW gluon distribution 
leads to much smaller cross section. The KL gluon distribution produces the 
hardest pt spectrum. Rather different slopes in transverse momentum of c 
(or c) are obtained for different UGDFs. This differences survive after con- 
voluting the inclusive quark/antiquark spectra with fragmentation functions. 
Thus, in principle, precise distribution in transverse momentum of charmed 
mesons should be useful to select a “correct” model of UGDF. A detailed 
comparison with the experimental data requires, however, a detailed know- 
eledge of fragmentation functions. 

The inclusive spectra are not the best observables to test UGDF [T]. Let 
us come now to correlations between charm quark and charm antiquark. 

The azimuthal angle correlation is the most popular observable in this 
context. There is a trivial relation between the relative azimuthal angle dis¬ 
tribution and the average transverse momenta of gluons. This is illustrated in 
Figinifor the simple Gaussian smearing with different values of the parameter 
(Jo ^ . In general, the larger (Jq the more decorrelation can be observed. In the 
literature both on-shell (see e.g.|3]) and off-shell (see e.g. irnsi) kinematics 
is used. Our calculation shows that there is almost no difference whether 
on-shell or off-shell kinematics is applied. In FigElwe show the angular cor¬ 
relations for the Kwiecihski UGDF for different values of the ho parameter. 
The effect here is smaller than for the simple Gaussian smearing. This is 
due to the fact that here the perturbative broadening is included explicitly 
and the parameter ho takes into account only nonperturbative part. The use 
of on-shell matrix elements instead of off-shell ones leads to slightly stronger 
back-to-back correlations. 

In FiglSlwe compare results for different unintegrated gluon distribution 
from the literature. Quite different results are obtained for different UGDF. 
The nonperturbative GBW glue leads to strong azimuthal correlations be¬ 
tween c and c. In contrast, BFKL dynamics leads to strong decorrelations 
of azimuthal angles of charm and anticharm quarks. The saturation-idea 
inspired KL distribution, as well as BFKL and KMR distributions lead to 

^The range in azimuthal angle is somewhat artificially increased to (-360°,360°) 
to better visualize the details of the distributions. The normalization is such that 
/- 27 r d(() = a. 
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an local enhancement for 0 c 5 ~ 0 which is probably due to gluon splitting s- 
channel subprocess. In the last case there is a sizeable difference between the 
result obtained with on-shell (left panel) and off-shell (right panel) matrix 
elements. Gluon distribution obtained by the KMR method and the KL and 
BFKL distributions generate very similar shapes in azimuthal angle. The 
gluon distribution obtained by the solution of the Kwiecihski equation gives 
somewhat narrower distribution around 0 = ivr. 

The initial state gluon transverse momentum is also responsible for the 
decorrelation of the charm quark, charm antiquark transverse momenta. In 
the case of absence of any gluon transverse momenta in leading-order collinear 
approximation pi * = p 2 ,t- The effect of transverse momenta is shown in FigEl 
for different unintegrated gluon distributions from the literature. Although 
the main strength is concentrated along the diagonal pi^t = P 2 ,t, there is a 
smearing over whole two-dimensional space of pi^t x P 2 ,t- As for the case of 
azimuthal correlations, different pattern can be observed for different UGDF. 

In general, the shape of the two-dimensional map is governed by the 
UGDF used and the matrix element squared. It is of interest to unfold these 
two effects. In FiglTUl we show an average matrix element squared (obtained 
with the help of BFKL UGDF) in the pi^t x P 2 ,t space. A rather weak de¬ 
pendence can be observed which means that UGDF is the crucial element 
responsible for the variation of the cross section in the pi^t x p 2 ,t space. We 
compare the on-shell (left panel) and off-shell (right panel) matrix elements 
squared averaged with UGDF over the phase space. In contrast to the more 
inclusive cases considered above, here the difference is quite sizeable. In 
Fig^^we present the corresponding ratio of both two-dimensional functions. 

What is the origin of populating the asymmetric hnal state conhgura- 
tions? The regions pi^t P 2 ,t and pi^t ^ P 2 ,t are populated by asymmetric 
initial conhgurations ^ -C t- This is demonstrated in FiglT^ 

where we show average value of kI ^ and t obtained for example with KMR 
gluon distributions. The smallest values of ^ and K 2 1 are obtained along 
the diagonal {pi^t = P 2 ,t)- This means that along diagonal the soft nonper- 
turbative part of UGDF is tested. The larger distance from the diagonal, 
the larger ^ and t are sampled. The KMR is known to have long tails in 
gluon transverse momenta. This in conjuction with Fig^J implies relatively 
large cross section for pi^t ^ P 2 ,t or pi^t -C P 2 ,t- 

In Fig El we show two-dimensional maps in pi^t and p 2 ,t obtained for four 
different complementary regions of ^ and 

(a) < 10 GeV2, ^2^ < 10 GeV^, 

(b) 4t < 10 GeV2, 4t > 10 GeV^, 

(c) > 10 GeV^, < 10 GeV2, 

(d) 4t > 10 GeV2, 4t > 10 GeVT 
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The asymmetric (in gluon transverse momentum) configurations contribute 
a large fraction to the integrated cross section. In Table 1 we ensambled 
fractions of different regions of k1 x kI- For completeness we placed in the 
table the integrated cross section for the cc pair production. The KMR and 
KL have largest fractions of asymmetric initial state configurations. The 
almost vanishing fraction of asymmetric configurations for GBW glue is due 
to a lack of pQCD effects. A rather schematic KL distribution includes 
the hard components only in a simple parametric form. Our purely model 
discussion here is a bit academic because initial state configurations are not 
directly the observables. Therefore the effect of the large tails in UGDFs can 
be only tested indirectly by studying asymmetric configurations in the final 
state (charmed mesons, muons) space pi^t x P 2 ,t- 

The decorrelation in pi ^ and p 2 ,t can be also studied in the collinear 
next-to-leading order approach. It is worth to stress here, however, that 
collinear approach even in higher order approximation is unreliable in the 
region pi ^ p 2 ,o In order to avoid conceptual and/or numerical problems 
asymmetric cuts on jets transverse momentum are usually assumed. On 
experimental side this means removing a big portion of statistics from the 
correlation studies. Our approach is free of the mentioned above problems 
and can be applied on the whole plane pi^t x p 2 ,t- 

In principle, there are many more correlation observables possible having 
in view quite reach phase space for two particle production. Here we shall 
concentrate on one more observable - namely the azimuthal angle correlation 
for different regions of quark-antiquark transverse momenta. In FigOl we 
show such correlations but now with extra conditions on quark-antiquark 
transverse momenta. Generally, the larger transverse momenta the more 
back-to-back correlation is observed. The details depend, however, on a 
particular model of UGDF. 

In FigElwe show similar correlations for very asymmetric values of pi^t 
and p 2 ,t- In this domain the cross section obtained with BFKL, KL and KMR 
UGDFs is much larger than the cross section obtained with the Kwiecihski 
UGDF. While the first three UGDFs lead to almost complete decorrelation of 
c and c, for the last UGDF we observe clear back-to-back correlations. Fur¬ 
thermore the absolute cross section in the latter cases is much smaller. This 
distinct difference causes that this observable is very promissing in testing 
the QGD dynamics. 

In order to shed more light on the relation between the c/c transverse 
momentum and azimuthal correlations in Figwe show the effect of trans¬ 
verse momentum of c or c and relative azimuthal angle on a two-dimensional 
map {pt, 4>-), where pt is transverse momentum of either quark or antiquark. 
The same effect as in Fig HU is shown now in a continuous way. There is also 
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a technical lesson from the inspection of FiglTHl The fignre shows that one 
needs somewhat more refine integration in for larger valnes of pi^t and 
P2,t- 

In the present exploratory calcnlation we have stndied correlations be¬ 
tween charm qnark and charm antiquark. In practice one measures either 
muons or charmed mesons. It seems that a study of x pt{D*~) or 

Pt{p~^) X pt{p~) or mixed pt{D*^) x pt{p~) or pt{p^) x pt{D*~) correlations 
would be optimal to test UGDFs. This will be a subject of our next stud¬ 
ies. It is not clear to us if such measurements are feasible with the present 
Tevatron apparatus. 

5 Conclusions 

We have compared quantitatively different methods to include gluon trans¬ 
verse momenta and their effect on inclusive spectra as well as on cc correla¬ 
tions. Different UGDFs from the literature were used. 

The gluon-gluon luminosity factor is more strongly varying factor over the 
phase space than the matrix element squared. Gonsequently, the shape of 
inclusive spectra only weakly depends on whether on-shell or off-shell matrix 
element is used. In contrast there is a sizeable effect for more exclusive 
variables, dependent e.g. on the position in the pi^t x P 2 ,t space. 

We have shown that the analysis of azimuthal correlations and decor¬ 
relations on the Pi,t{c transverse momentum) x P 2 ,t{c transverse momen¬ 
tum) plane as well as the combination of both are very promissing to verify 
UGDF. We have focussed on the region of pi^t 3> P 2 ,t and pi^t ^ P 2 ,t- There 
max{fi:i K 2 ,t} ~ niax{pi ^,^ 2 , 4 } and off-shell matrix element must be used. 
In this region the shape of azimuthal correlation function depends on both 
UGDF and matrix element. 

Although we have calculated correlation observables for charm quarks 
and antiquarks it is obvious that the main effects should be very similar 
for charmed mesons and/or muons. The correlations between D*^{cd) and 
D*~ (cd) or /i+ and are the most promissing in this respect as they reflect 
relatively well the cc final state. Estimating the rates of events with one 
D meson with small transverse momentum and second D meson with large 
transverse momentum would be a useful test of the QGD dynamics. The 
analysis of meson correlations will be a subject of our next analysis. 
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Table 1: Fractions in % of contributions from different regions of x fo 
the total cross section for different UGDFs from the literature for W = 1.96 
TeV, with off-shell kinematics and off-shell matrix elements. In the case of 
Kwiecihski UGDF: bo = 1 GeV“^, fx'jp = A in the case of KMR UGDF: 
= A ml. Below for brevity = 10 GeV^ is assumed. 


UGDF 


ws^sKs^m 



ws^sKs^m 

Kwiecihski 

0.85 

91.87 

3.99 

3.99 

0.15 

KMR 

0.55 

52.43 

19.83 

19.83 

7.91 

BFKL 

1.44 

59.65 

17.26 

17.26 

5.83 

GBW 

0.12 

100.0 

0.001 

0.001 

0.0 

KL 

0.51 

65.49 

15.37 

15.37 

3.76 


TABLES 
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Figure 1: Inclusive da/dpt for charm/anticharm production at W = 1.96 TeV 
in the approach of naive Gaussian smearing of gluon transverse momenta. 
The calculation with different values of Uo = 0.25 (dotted), 0.50 (dashed), 
1.0 (solid), 2.0 (dash-dotted) GeV are shown. The results with on-shell 
kinematics are shown in panel (a) and results with off-shell kinematics in 
panel (b). In this calculation both factorization and renormalization scales 
were hxed for 4 
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Figure 2: Inclusive da/dpt for charm/anticharm production at W = 1.96 
TeV for the Kwiecihski unintegrated gluon distributions for = 0.5, 1.0, 2.0 
GeV“^ and for = 4m^. The results with on-shell matrix element is shown 
in the left panel and results with off-shell matrix element in the right panel. 
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Figure 3: Inclusive da/dpt for charm/anticharm production at W = 1.96 
TeV for the Kwiecihski UGDF for different factorization scales 
(dashed), 4 ml (solid). In this calculation ho = 1.0 GeV“^. 
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Figure 4: Inclusive da/dpt for charm/anticharm production at W = 1.96 TeV 
for the Kwiecihski UGDF for on-shell (dashed) and off-shell (solid) matrix 
element. In this calculation 6o = 1-0 GeV“^, /i^ = = 4 m^. 
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Figure 5: Inclusive da/dpt for charm/anticharm production at W = 1.96 
TeV for different UGDFs. The meaning of the curves is as follows: solid - 
Kwiecihski, thick dashed - KL, thin dashed - GBW, dash-dotted - BFKL, 
dotted - KMR. 




Figure 6: Azimuthal angle correlations for different values of the Gaussian 
parameter uo = 0.25, 0.50, 1.0, 2.0 GeV for on-shell (panel a) and off-shell 
(panel b) kinematics. Here 0 < Pi,t,P 2 ,t < 10 GeV. 
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Figure 7: Azimuthal angle correlations for the Kwiecihski UGDF for different 
values of the parameter bo = 0.5, 1.0, 2.0 GeV“F In this calculation /i|, = 
= 4 was taken. 
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Figure 8: Azimuthal angle correlations for different UGDF in the literature: 
Kwiecihski (solid), BFKL (dash-dotted), GBW (thin dashed), KL (thick 
dashed) and KMR (dotted). 
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Figure 9: Two-dimensional distributions in pi ^ of charm quark and p 2 ,t of 
charm antiquark for different UGDF in the literature: (a) KMR (/i|. = 4m^), 
(b) Kwiecihski {bo = 1.0 GeV“^, f^p = 4m^), (c) BFKL, (d) KL. 
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Figure 10: Average matrix element squared as a function of of charm 
quark and p 2 ,t of charm antiquark for on-shell (left panel) and off-shell (right 
panel) matrix element. In this calculation BFKL UGDF was used. 
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Figure 12: Average transverse momenta of initial gluons or as a 
function of pi^t and p 2 ,t- In this calculation Kimber-Martin-Ryskin UGDF 
was used. 
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Figure 13: Two-dimensional distributions in pi^t of charm quark and p2,t 
of charm antiquark for the KMR UGDF and for different regions of gluon 
transverse momenta specihed in the text. 
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Figure 14: da for different intervals of charm and anticharm transverse 
momenta: 

0.0 GeV < < 2.5 GeV (solid), 

2.5 GeV < pi^t,P 2 ,t < 5.0 GeV (dashed), 

5.0 GeV < pi,t,P 2 ,t < 7.5 GeV (dash-dotted), 

7.5 GeV < Pi,t,P 2 ,t < 10.0 GeV (dotted) 

for different UGDF: (a) KMR (/i^ = 4m^), (b) Kwiecihski {bo = 0.5 GeV“\ 
^2 = 4^2)^ BFKL, (d) KL. 
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Figure 15: da/d(f)_ for 

0.0 GeV < Pi t < 2.5 GeV and 7.5 GeV < p 2 ,t < 10. GeV or 
7.5 GeV < pi^t < 10. GeV and 0.0 GeV < p 2 ,t < 2.5 GeV 
for different UGDF: KMR = 4m^) - dotted, Kwiecinski {bo 
p?p = 4ml) ~~ solid, BFKL - dash-dotted, and KL - dashed. 
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